Virial Theorem

In mechanics, the virial theorem provides a general equation that relates the average over time of the total kinetic energy of a stable system of discrete particles, bound by potential forces, with that of the total potential energy of the system. Mathematically, the theorem states
$\left\langle T \right\rangle = -\frac12\,\sum_^N \bigl\langle \mathbf_k \cdot \mathbf_k \bigr\rangle$
where is the total kinetic energy of the particles, represents the force on the th particle, which is located at position , and angle brackets represent the average over time of the enclosed quantity. The word virial for the right-hand side of the equation derives from ''vis'', the Latin word for "force" or "energy", and was given its technical definition by Rudolf Clausius in 1870.

The significance of the virial theorem is that it allows the average total kinetic energy to be calculated even for very complicated systems that defy an exact solution, such as those considered in statistical mechanics; this average total kinetic energy is related to the temperature of the system by the equipartition theorem. However, the virial theorem does not depend on the notion of temperature and holds even for systems that are not in thermal equilibrium. The virial theorem has been generalized in various ways, most notably to a tensor form.

If the force between any two particles of the system results from a potential energy that is proportional to some power of the interparticle distance , the virial theorem takes the simple form
$2 \langle T \rangle = n \langle V_\text \rangle.$

Thus, twice the average total kinetic energy equals times the average total potential energy . Whereas represents the potential energy between two particles of distance , represents the total potential energy of the system, i.e., the sum of the potential energy over all pairs of particles in the system. A common example of such a system is a star held together by its own gravity, where equals −1.

History

In 1870, Rudolf Clausius delivered the lecture "On a Mechanical Theorem Applicable to Heat" to the Association for Natural and Medical Sciences of the Lower Rhine, following a 20-year study of thermodynamics. The lecture stated that the mean vis viva of the system is equal to its virial, or that the average kinetic energy is equal to the average potential energy. The virial theorem can be obtained directly from Lagrange's identity as applied in classical gravitational dynamics, the original form of which was included in Lagrange's "Essay on the Problem of Three Bodies" published in 1772. Karl Jacobi's generalization of the identity to ''N'' bodies and to the present form of Laplace's identity closely resembles the classical virial theorem. However, the interpretations leading to the development of the equations were very different, since at the time of development, statistical dynamics had not yet unified the separate studies of thermodynamics and classical dynamics. The theorem was later utilized, popularized, generalized and further developed by James Clerk Maxwell, Lord Rayleigh, Henri Poincaré, Subrahmanyan Chandrasekhar, Enrico Fermi, Paul Ledoux, Richard Bader and Eugene Parker. Fritz Zwicky was the first to use the virial theorem to deduce the existence of unseen matter, which is now called dark matter. Richard Bader showed the charge distribution of a total system can be partitioned into its kinetic and potential energies that obey the virial theorem. As another example of its many applications, the virial theorem has been used to derive the Chandrasekhar limit for the stability of white dwarf stars.

Illustrative special case

Consider particles with equal mass , acted upon by mutually attractive forces. Suppose the particles are at diametrically opposite points of a circular orbit with radius . The velocities are and , which are normal to forces and . The respective magnitudes are fixed at and . The average kinetic energy of the system is

$\langle T \rangle = \sum_^N \frac12 m_k \left, \mathbf_k \^2 = \frac12 m, \mathbf_1, ^2 + \frac12 m, \mathbf_2, ^2 = mv^2.$
Taking center of mass as the origin, the particles have positions and with fixed magnitude . The attractive forces act in opposite directions as positions, so . Applying the centripetal force formula results in:
$-\frac12 \sum_^N \bigl\langle \mathbf_k \cdot \mathbf_k \bigr\rangle = -\frac12(-Fr - Fr) = Fr = \frac \cdot r = mv^2 = \langle T \rangle,$
as required. Note: If the origin is displaced then we'd obtain the same result. This is because the dot product of the displacement with equal and opposite forces , results in net cancellation.

Statement and derivation

Although the virial theorem depends on averaging the total kinetic and potential energies, the presentation here postpones the averaging to the last step.

For a collection of point particles, the scalar moment of inertia about the origin is defined by the equation
$I = \sum_^N m_k \left, \mathbf_k \^2 = \sum_^N m_k r_k^2$
where and represent the mass and position of the th particle. is the position vector magnitude. The scalar is defined by the equation
$G = \sum_^N \mathbf_k \cdot \mathbf_k$
where is the momentum vector of the th particle. Assuming that the masses are constant, is one-half the time derivative of this moment of inertia
$\begin \frac12 \frac &= \frac12 \frac \sum_^N m_k \mathbf_k \cdot \mathbf_k \\ &= \sum_^N m_k \, \frac \cdot \mathbf_k \\ &= \sum_^N \mathbf_k \cdot \mathbf_k = G\,. \end$
In turn, the time derivative of can be written
$\begin \frac & = \sum_^N \mathbf_k \cdot \frac + \sum_^N \frac \cdot \mathbf_k \\ & = \sum_^N m_k \frac \cdot \frac + \sum_^N \mathbf_k \cdot \mathbf_k \\ & = 2 T + \sum_^N \mathbf_k \cdot \mathbf_k \end$
where is the mass of the th particle, is the net force on that particle, and is the total kinetic energy of the system according to the velocity of each particle
$T = \frac12 \sum_^N m_k v_k^2 = \frac12 \sum_^N m_k \frac \cdot \frac.$

Connection with the potential energy between particles

The total force on particle is the sum of all the forces from the other particles in the system
$\mathbf_k = \sum_^N \mathbf_$
where is the force applied by particle on particle . Hence, the virial can be written
$-\frac12\,\sum_^N \mathbf_k \cdot \mathbf_k = -\frac12\,\sum_^N \sum_^N \mathbf_ \cdot \mathbf_k \,.$

Since no particle acts on itself (i.e., for ), we split the sum in terms below and above this diagonal and we add them together in pairs:
$\begin \sum_^N \mathbf_k \cdot \mathbf_k & = \sum_^N \sum_^N \mathbf_ \cdot \mathbf_k = \sum_^N \sum_^ \left( \mathbf_ \cdot \mathbf_k + \mathbf_ \cdot \mathbf_j \right) \\ & = \sum_^N \sum_^ \left( \mathbf_ \cdot \mathbf_k - \mathbf_ \cdot \mathbf_j \right) = \sum_^N \sum_^ \mathbf_ \cdot \left( \mathbf_k - \mathbf_j \right) \end$
where we have assumed that Newton's third law of motion holds, i.e., (equal and opposite reaction).

It often happens that the forces can be derived from a potential energy that is a function only of the distance between the point particles and . Since the force is the negative gradient of the potential energy, we have in this case

$\mathbf_ = -\nabla_ V_ = - \frac \left( \frac \right),$

which is equal and opposite to , the force applied by particle on particle , as may be confirmed by explicit calculation. Hence,

$\begin \sum_^N \mathbf_k \cdot \mathbf_k &=\sum_^N \sum_^ \mathbf_ \cdot \left( \mathbf_k - \mathbf_j \right) \\ &=-\sum_^N \sum_^ \frac \frac \\ &=-\sum_^N \sum_^ \frac r_. \end$

Thus, we have
$\frac = 2 T + \sum_^N \mathbf_k \cdot \mathbf_k = 2 T - \sum_^N \sum_^ \frac r_.$

Special case of power-law forces

In a common special case, the potential energy between two particles is proportional to a power of their distance
$V_ = \alpha r_^n,$
where the coefficient and the exponent are constants. In such cases, the virial is given by the equation
$\begin -\frac12\,\sum_^N \mathbf_k \cdot \mathbf_k &=\frac12\,\sum_^N \sum_ \frac r_ \\ &=\frac12\,\sum_^N \sum_ n \alpha r_^ r_ \\ &=\frac12\,\sum_^N \sum_ n V_ = \frac\, V_\text \end$
where is the total potential energy of the system
$V_\text = \sum_^N \sum_ V_ \,.$

Thus, we have
$\frac = 2 T + \sum_^N \mathbf_k \cdot \mathbf_k = 2 T - n V_\text \,.$

For gravitating systems the exponent equals −1, giving Lagrange's identity
$\frac = \frac12 \frac = 2 T + V_\text$
which was derived by Joseph-Louis Lagrange

HOME


Content is Copyleft
Website design, code, and AI is Copyrighted (c) 2014-2017 by Stephen Payne


Consider donating to Wikimedia



As an Amazon Associate I earn from qualifying purchases